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Abstract. We consider diffraction at random point scatterers on generai discrete point 
sets in R"^, restricted to a finite volume. We allow for random amplitudes and random 
dislocations of the scatterers. We investigate the speed of convergence of the random 
scattering measures applied to an observable towards its mean, when the finite volume tends 
to infinity. We give an explicit universal large deviation upper bound that is exponential in 
the number of scatterers. The rate is given in terms of a universal function that depends on 
the point set only through the minimal distance between points, and on the observable only 
through a suitable Sobolev-norm. Our proof uses a cluster expansion and also provides a 
central limit theorem. 
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1. INTRODUCTION AND SETUP 

The study of the diffraction theory of 'ordered point sets' is a classical subject to 
physicists: Crystals produce sharp diffraction images, with bright spots known as Bragg 
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peaks. It is known since the eighties (and mathematicaUy weU-understood by now) 
that this is also true for quasi-crystals ([Hof95b]). They posses fong-range order but no 
translation symmetry and show geometrically intriguing diffraction patterns. Of course, 
sharp behavior of the scattering image occurs only in the limit of an infinite system of 
scatterers, much in analogy to the sharpness of phase transitions in statistical mechanics. 

What changes, if one is adding randomness ('disorder') to the picture? It is reason- 
able to expect that there should be a well-defined limit of the diffraction image when 
the number of scatterers tends to infinity, under natural assumptions. What assump- 
tions exactly do we need mathematically? Do sample fluctuations matter? Do we have 
control over corrections to the infinite volume behavior when the number of scatterers 
is finite? There have been few mathematical papers about the first two questions (see 
however [BaaMoo98] , [BaaHoeOO], [Hof95a]), and no results at all about the finite volume 
behavior. 

Suppose at first one is choosing the scatterers according to some translation-ergodic 
distribution while keeping the position fixed on a perfect crystal. Then the scattering 
images will converge to their disorder-averages in a distributional sense, by soft ergodicity 
arguments. This is true for almost any realization of the scatterers, using the ergodic 
theorem. Note however that these arguments do not provide any control over the finite 
volume corrections of the observed scattering image. 

When one gives up translation invariance of the underlying structure or the distribu- 
tion, there are not very many mathematical results in the literature. Ergodicity argu- 
ments are not available any more and one must resort to explicit methods. Subjecting 
the sites of a quasicrystal to an i.i.d. thermal motion leads to an infinite volume pic- 
ture that is well-known from crystals: The intensity of the sharp peaks is reduced by a 
Debye- Waller factor with a diffuse background appearing. It is however not difficult to 
rigorously justify this kind of law-of-large number result when no control over the speed of 
convergence is required (see [Hof95a]). For mathematical results about the scattering at 
random tilings in the infinite volume limit we refer the reader to the recent review article 
of [BaaHoeOO]. This paper provides a number of interesting and pedagogical examples 
and can also serve as a good introduction to mathematical scattering theory. 

In the present paper we provide a contribution to the diffraction theory of random 
scatterers on general point sets by answering the question 

How large is the probability for a deviation of the scattering image of a finite 
portion of scatterers from its sample average? 

This is the probabilistically natural question for self-averaging and the experimental- 
ist's question 'Is my system large enough?'. Our emphasis in this paper is that we do not 
assume a lattice structure, quasiperiodic structure or any symmetry of the set or of the 
distribution. For our results we only assume a minimal distance between the points of 
the reference point set. In particular all of our results hold for lattices or quasi-crystals. 
The answer then immediately leads to results about convergence of the scattering images 
for almost every realization, when it is combined with information about the behavior of 
the mean (without unnecessary assumptions as found in the literature). 
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We consider two types of randomness: A) We choose the scattering amphtudes ac- 
cording to a random distribution while keeping their locations fixed or B) subject them 
to thermal motion around their sites. We can treat both types in a unified way. In order 
to investigate the selfaveraging properties in a physically meaningful way, we adopt the 
following point of view. Fix an observable, modelling the counter used in an experiment, 
and look at the result of the corresponding measurement. Then estimate the proba- 
bility distribution of the resulting quantity (that is random w.r.t sample fluctuations) 
in any finite volume. This provides much more information than the mere statement 
of the convergence of the scattering measures in the sense of distributions, for almost 
every realization of the scatterers, as it is conventionally done ([BaaMoo98], [BaaHoeOO], 
[Hof95a]). 

We note that, for a typical observable (e.g. Gaussian test function) the resulting ex- 
pression will involve all autocorrelation coefficients of the array of scatterers and there- 
fore does not trivially decompose into independent parts, even for scatterers that behave 
independently. In the language of statistical mechanics, the observable produces an inter- 
acting system! Treating the autocorrelation coefficients as individual random variables 
without using their dependence would suflBce for a mere convergence result, but would 
lead to very bad large-deviation estimates. To deal with this interaction it will there- 
fore be appropriate to employ (high-temperature) expansion methods from statistical 
mechanics, as will become clear soon. This shows the usefulness of such methods to 
give sharp explicit results, even in situations that do not a priori smell like dependent 
spin systems and Gibbs-measures. So, it would be nice if the paper could also serve as 
a motivation for probabilists and mathematical physicists who are sceptical about the 
use of expansions to take a closer look. In fact we restrict ourselves to the situation of 
independent scatterers to keep the technicalities down. To generalize the method to the 
case of weakly coupled scatterers is possible, but it would complicate the theorems, and 
make the general idea less transparent. We therefore leave it to a future paper. 

In this setup we will provide general upper bounds on the probability that the mea- 
surement in finite volume deviates from its mean, and even provide explicit numerical 
values of the constants appearing. These estimates are universal in the sense that they 
depend only: 1) on the minimal distance between sites, but not on details of the point set; 
and 2) on the concentration of the observable, measured in a suitable Sobolev-norm of its 
Fourier-transform, but not on any more details. The fact that the estimate depends on 
the point set only through the minimal distance is important because one might want to 
be able to interpret the diffraction images without knowing beforehand the geometrical 
structure of the point set, while having some physical a priori-estimate on the mimimal 
distance. 

Computing the average of a scattering image is simple and seeing whether it con- 
verges or not reduces to the knowledge of the autocorrelation structure already needed 
to understand the deterministic image of the point-set (see Appendix A). 

Let us now define the models and state our results more precisely. We discuss the 
scattering image at infinity that is created by single-scattering at (a finite collection of) 
the point-scatterers described by the following random measures. 
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Model A: Disordered scattering amplitudes, fixed positions 



We look at the complex random measure ( 'random Dirac comb ') given by 



X 



(1.1) 



x&r 



where 5x denotes the Dirac-measure at the site x. The point set FcM'^ is assumed to 
be countable. Here rjx are random, possibly complex scattering amplitudes, independent 
over the sites x eT. They are assumed to be bounded. 

Model B: Thermally dislocated scatterers 

Here we look at the scattering image of the random measure given by 



where Ux are random ('thermal') dislocations taking values in M^, independent over the 
sites X G r. They are assumed to be bounded, too. 

Fix any finite volume Fj-CF. Then, the object that contains all information about 
the scattering image of the points in F^. is the finite volume scattering measure which is 
the Fourier-transform of the corresponding finite volume autocorrelation measure. (For 
a summary of the basic notions of mathematical scattering theory, see Appendix A and, 
for more details, e.g. Chapter II of [BaaHoeOO].) Here, for Model A the autocorrelation 
measure in the finite volume Fy. is given by 



where the star denotes complex conjugate and the y-sum is over all difference vectors in 
F^. Since we allow F^. to be any finite set, we have chosen the natural normalization by 
the number of points (in contrast to [BaaHoeOO]). This leads to simpler formulas in our 
theorems. For Model B we put 



for the finite volume autocorrelation measure. Suppose now a measurement on the scat- 
tered intensity is performed that is described by an observable (p{k) in Fourier-space, 
modelling the coimter. Usually it is assumed to be a real Schwartz function. The corre- 
sponding result of the measurement is then given by 7^(y') = / ^r{k)(p{k)dk. Here the 
Fourier-transform of a tempered distribution 7 is defined by duality, 7(v) = 7(<^), where 
(f denotes the Fourier-integral of the Schwartz- function (f over W . (For a quick reminder 




(1.2) 



a;€r 




(1.3) 




(1.4) 
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of the explanation for this and some comments, see Appendix A. For more expository- 
details, see [BaaHoeOO].) The sample average of the measurement is J ii{dini)^'!^{ip). This 
object has the correct normalization (by the total number of scatterers) to be able to 
converge to a well-defined limit, as it does of course for non-random scatterers on a 
lattice. (We remark that this average won't converge [e.g. for i.i.d. scatterers] along 
any sequence of volumes F^, but e.g. an increasing sequence of balls or cubes works for 
crystals and quasi-crystals.) 



Main result for Model A 

Let us formulate the bound in the simplest form, which is suitable for the computation 
of explicit numbers bounding the probability of a large deviation. It makes explicit 
the uniformity of the large deviation upper-bound, independently of the set F (other 
than the minimal distance between points in F), the form of the distribution (other 
than through uniform bounds on the magnitude of the scatterers), and the observable ip 
(other than through a Sobolev-norm) . We also provide a different version of the large 
deviation estimate in Chapter 4 under the name 'Addition to Theorem 1 '. It is slightly 
sharper in certain cases but less useful for direct application. In Chapter 4 we also give 
a corresponding Central Limit Theorem. 

Now, to state the theorem we define the following Sobolev-norm involving integrals of 
derivatives up to the order of the dimension i/, where we also introduce a scaling factor 
a/ 2. For a function ^ : — > C we put 

Here |i?i| denotes the volume of the z/-dimcnsional unit ball. The symbol 
d'^g{y) : (R'^)'^ — >■ R'^ denotes the k-th differential of g at the point y and 
ll'^'^fi'd/) II '■— sup|^,^|= \d''g{y)[vi, . . . ,Vk]\ is the usual norm of a /c-multilinear map- 

ping, at any fixed point y, where |f | denotes the Euclidean norm. 
Then we have 

Theorem 1. Suppose that F^cR"^ is any finite set and denote the minimal distance be- 
tween its points by a. Assume that rj = {r]x)x€r Q'^e (possibly complex) random variables, 

independent, but not necessarily identically distributed. Denote their distribution by /i. 
Suppose the uniform bounds l/J-i^x)] < M < oo and \r]x — A*(%)l ^ B < oo, for all 
X e Fr, for jjL-a.e. realization. 

Then the corresponding random scattering image 7^(9?) in the finite volume F^ obeys 
the universal large deviation estimate 

l^ ( 7?(^) - / /^(^^)7."(^) >e^< 2exp(-|F,| x (1.6) 
for any £ > 0, for any function (/? : R'^ 1— > R, s.t. its Fourier-transform has finite norm 

\\'Ph,a- 
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Here K = 2MB + B^, and J is a nonnegative, convex, strictly monotone function 
that is independent of the form of the distribution /i and the set T. 

Remark. The universal function J : [0, oo) — > [0, cxo) in the theorem has the form 

\ d{e-d{2 + dD)) , else 
where the numerical constants can be chosen like d 0.0525, and D = 4.54 • 10^. Note 

-2 

the asymptotics J{e) ~ ^ for £ ], and J{e) ~ de for e f oo. 

A structural and practical virtue of the form of the large-deviation estimate of Theo- 
rem 1 lies in the fact that its dependence on the observable ip is formulated entirely in 
terms of the continuum-object ||<^||i/,a- All details of the set have disappeared! This 
Sobolev-norm can be computed (at least numerically) with little effort, and so one may 
easily derive explicit numbers. 



Remark. Note the natural fact that the bound is scale-invariant in the following way: 
Suppose the counter is modelled by a probability density tpaik) = (pi{k / a) in Fourier- 
space with variance ('precision of measurement') cr^. (Think e.g. of a Gaussian!) Then, 
by scaling we have ||(^a|U,a = ll'^ilU.aa- So we have ||<^a|U,a ~ (ao-)"''/jg, \'pi{y)\dy with 
(T I 0, when a is fixed (under the condition that the higher derivatives arc intcgrable). 
This immediately controls the deterioration of our large deviation estimate when we make 
a smaller to increase the precision of measurement of the scattering image. (Without loss 
we could have chosen our length-scale in such a way that a = 1 from the beginning, so 
that the general statement is regained by rescaling the observable in /c-space. We believe 
however that the present form of the theorem is more intuitive.) 

Remark. The norm appearing is finite in particular for the commonly used Schwartz- 
test-functions. So our result in particular implies convergence-statements in the sense 
of tempered distributions. Suppose we are given an increasing sequence of volumes F^. 
Then we immediately obtain the strong law of large numbers as a consequence of Theorem 
1, saying that the centered autocorrelation measure applied to a test function (p whose 
Fourier transform has finite norm converges to zero, for P-a.e. i]. This follows trivially by 
summing the exponential bound (1.6) over the volumes using the Borel-Cantelli Lemma. 

Remark. The fact that the dependence on the observable on the point-set F and on 
the distribution can be expressed in terms of the handy quantity is not a 

priori obvious. The occurrence of the norm however is not difficult to understand. It can 
be motivated by noting that 4(i^ ||(^||i.^a)^/|Fr | is an upper bound for the ^u-expectation 
of the square of the modulus inside the probability on the l.h.s. of (1.6). (This is seen 
using the independence of the scatterers by Fourier-transform, and substituting the norm- 
estimate of Proposition 3.) Believing in the corresponding Gaussian behavior, the small 
£-behavior given in (1.7) should follow. An essential part of the real proof consists in 
estimating all the higher moments contained in the Laplace-transform in terms of powers 
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of Kllc^lliy^a- One can not expect a large deviation principle [which would in particular 
mean the existence of the limit — lim^ log// (| ... | > s)] without any assumptions on 
the set r other than minimal distance. In fact, without further assumptions on F, the 
Laplace-transform won't converge. 



Main result for Model B 

For the Model B of thermal dislocations the result is quite analogous. Here, however 
the Sobolev-norm of the variation of the Fourier-transform of the observable appears. 
Again, there will be a sharper version of this result in Chapter 4 that is called 'Addition 
to Theorem 2', and the Central Limit Theorem. 

Theorem 2. Suppose again that TrdW is any finite set and denote the minimial dis- 
tance between its points by a. 

Suppose that the dislocations uj = {ujj:)nc^rr ho,ve independent, not necessarily identical 
distribution n, such that \u)x\ < 5 < a/4, for all x G F^, fi-a.s. 

Then the finite volume scattering image 7r (v) obeys the universal large deviation 
estimate 

,i\^n^)-j,(^m,f)\>e) <2exp(-|r.|xj(j^^^^— -)) (1.8) 

The function J has the same form as the function J from Theorem 1 (see (i-7)), but 
with the slightly better constant D = 4.38- 10^ (< D) instead of D, and the same constant 
d. 

Here the appearing norm has the obvious meaning obtained by extending the previous 
definition (1.5) that was given for functions to linear functionals. It equals 
\\d9\\u,a = ^ ELo M ia/2y-'' /r- W^^^^ MWdy ■ 



Remark. The restriction S < a/4 is only for simplicity. The more general statement of 
the Addition to Theorem 2 stays true for any finite 6. The fact that the estimate involves 
the Sobolev-norm of the derivative of (p rather than the Sobolev-norm of (p itself is due 
to the fact that 7^(0;) is non-random when the function a is a constant. 

Remark. Note again the scale-invariance of the estimate, where of course the spatial 
distance d must be rescaled, too: As for Model A, take a rescaled observable (fa{k) = 
a~'^(^i{k/a) in Fourier-space. Then we have d ||(i<^(j||;/,a-4(5 = Sa ||(i(^i||;/,a(j-4(5o-- So, 
assuming that the higher derivatives are integrable, the quantity appearing in the large 
deviation estimate behaves like 

S \\d'fia\\u,a-46 ~ S{a - 45)~''cr~^'^+^) J^, \\d(pi{y)\\dy with cr J, and a, S fixed. 
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We conclude this introduction with an outhne of the rest of the paper along with 
some ideas of the proof. In Chapter 2 we derive bounds on the Laplace transform of 
the centered random scattering measures, applied to some observable ip. To do so we 
look at this quantity as a (possibly complex) Hamiltonian of a spin-system. Here the 
random variables modelling the scatterers (resp. their dislocations) play the role of 
spins. The Laplace-transform then becomes a partition function. We can treat it by 
high-temperature expansion methods from statistical mechanics, under the assumption 
that the interaction be small. The smallness of the interaction of the spin-system we 
need for the expansion will be guaranteed by smallness of the Fourier-transform in a 
suitable norm. From the point of view of the expansion it is natural to introduce discrete 
T-dependent norms, so that we can control the terms of the expansion with constants that 
are independent of the structure of F. The resulting bounds for the Laplace-transforms 
including the computation of numerical constants are provided in Proposition 1 for Model 
A. In Chapter 3 the work of Chapter 2 is adapted to treat Model B. In Chapter 4 we 
state the sharpened results of the 'Additions to Theorems 1 and 2' and the Central Limit 
Theorem, along with their proofs. They follow immediately from the norm-estimates on 
the Laplace-transform. 

In Appendix A we recall the basic notions of scattering theory for point scatterers. In 
Appendix B we give estimates on our discrete F-dependent norms in terms of Sobolev- 
norms that depend on F only through the minimal distance. This input is needed to 
show the uniformity in F and the nice bounds given in Theorems 1 and 2. 



2. NORM BOUNDS ON MOMENT GENERATING FUNCTION 

In this chapter we use an expansion to derive bounds on the Laplace transform of the 
random variable in question. We will look at this random variable as a Hamiltonian of a 
spin system. We will formulate the bounds obtained in this chapter in terms of a suitable 
discrete norm that is close to what is needed for the proof of convergence, and compute 
numerical constants. These constants are obtained employing the known Kotecky-Preiss 
estimates for abstract polymer models. The result of this is found in Proposition 1. 

Now, let us use the short notation 



for the nonnormalized centered autocorrelation measure applied to the function a '.W^ ^ 
C. This is the random variable in question. To derive bounds on its Laplace-transform 
and control the terms higher than second order in a we need a suitable norm. It turns 
out that the appropriate norm is the discrete Z^-type norm 




(2.1) 




(2.2) 
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Of course, for F = aU^ this is just an Z-'^-norm, for general F it is slightly more complicated. 
The result of this chapter is then the following. 

Proposition 1. Suppose that the scatterers rjx have independent, not necessarily identi- 
cal distribution /x, such that |// {rjx) | < M < cxo and \r]x — A*(??a;)| < B < oo for all x eT, 
jjL-a.s. 

Then there are universal constants d > ^, and D < oo, independent of the set F 

and the distribution (for given K), such that, whenever a : C is such that 

\o({x)\ = |q;(— for all a; e F, and i^||Q;||r < d, we have the estimate 



log/.(e^'^(«))-l/.(X,(a)f 



< |F^|L> X (i^||a||r)^ (2.3) 



with K := 2MB + B^. 

The values of the constants can be chosen like d — 0.0525, and D = 4.54 • 10^. 

Remark. Note that the quadratic term in oi under the modulus may not have a limit with 
r t oo, for general sets F, even in the i.i.d. case. Much less need the higher moments of 
Xr{a) possess a limit. The essential point is however that all higher order terms in a 
are estimated uniformly in the set F. This uniformity follows from the cluster expansion 
error bounds and some explicit work. 

Proof. We interpret X^{a) as the (negative) Hamiltonian of a spin-system with spin- 
variables Tjj:, X eVj. and open boundary conditions. It is then most intuitive from the 
point of view of statistical mechanics to write it in the form 



Xr{a) = Yl U^,^ (2 4) 



with the single-site potential 

Kc = a(0)(|?7a.p-/^(|%|')) (2.5) 

and the pair potential 

Ux,y =a{x - y) (r]xr]*y - f^iVxV^)) + «(y - x) (r]ynl - fi{VyV*x)) (2-6) 

for X y. Note that when a is the Fourier-transform of a real function, the pair 
potential is real. In general we allow it to be complex. Note that the potentials V and 
U are linear in the function a. Note that in general the interaction will act between 
all pairs scatterers 1]^. We look at the logarithm of the Laplace transform of Xr{a) 
which becomes the partition function of the spin-system. Then we want to compute the 
partition function of corresponding spin model to quadratic order in the strength of the 
interaction and control the remainder term. If we restrict ourselves to sufficiently 'small' 
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a this will allow us to perform a cluster expansion, as we will see. This corresponds to 
'small inverse temperature'. Such an expansion is in principle well-known in statistical 
mechanics, but we have to be careful about the precise assumptions we need on a and 
keep track of the constants appearing. It turns out that all we need to for the control is 
the quantity i^||Q:||r- 

Let us start. We note that 

IIK.II00 < |«(0)|S2, WU^^yWoo < 2K\a{x - y)\ (2.8) 

Then we put A*[V^](-) := A*( • e^^e^^ ^"')///(e^^sr^ to separate the single-site contribu- 
tions and write 



\{a;,2}cr^ / 



(2.9) 



EH n 

TcBr \{x,z}eT 



The set Br describes the set of edges on the complete graph with vertices F^. We write 
T = PiU ■ ■ - U Pn for the unique decomposition into connected components and call the 
Pi's polym,ers. A polymer P is thus of the form P = {{xi, ^i}, {x2, -22}, • • • , {a^fc, -Zfc}} 
and will be considered as a connected graph. There is the obvious notion of pairwise 
compatibility: Pi,P2 are compatible iff they don't have any sites in common. So we 
write the last expression as a sum over pairwise compatible families of polymers with 
polymer-activities that depend on a. 

n 

^ (e^^(«)) =11^ (e^^) X Yl Up^^i^) (2.10) 
xerr- (Pi,...,p„),i=i 

Here the polymer activity of a polymer is given by 

PP ^ pp{a) = ^^[V] I n (e^--l) I (2.11) 




This is the general formulation of a polymer partition function in an abstract polymer 
model. We want to perform the corresponding cluster-expansion for the logarithm of 
it. This is nothing but the Taylor-expansion when the polymer-activities are treated as 
independent (complex) variables pp. 

After this is done, we expand the activities pp to quadratic order as functions of 
a. Expanding the exponential in powers of a gives the following. Let us write gi{s) = 
^iLi jr the remainder term of the Taylor-series of the exponential and use that 
< 9li\s\). 
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For a general polymer we have the bound 

\pp{a)\< Yl gi{u^,z) (2.12) 
{x,z}eP 

with the abbreviation 

UxjZ '•— ||^a;,a:||oo (2.13) 

using the uniform bound on the scatterers. 

This is the bound we use to prove convergence of the expansion. We apply (a slight 
extension of) Proposition A.l that can be found in [KOI]. It says 

Proposition A.l of [KOI]. Suppose that ^(p^ p ) nr=i /^^i ^ polymer partition 
function, where: 'Polymers' P are graphs on a set Vj. having at least one edge. Two 
polymers are called compatible if they have disjoint vertex sets. The sum is over pairwise 
compatible families of polymers. Assume that the (possibly complex) activities pp satisfy 
the bounds 

\pp\ < e'^bep-'b where A := sup e""^"-^ < A* f=:i 0.110909 ^2 14) 

for some function Tb — Tx,y > on the set of edges on Fr, where the above b-sum is over 
all edges of the graph P. 

Then, the cluster expansion converges, i.e. the Taylor-series of the logarithm of the 
partition function has the representation 

n 

IlpPi = J2^(^ (2-15) 

(Pl,...,P„)ei=l C 

where the sum is over indecomposable subsets CcV . 'Indecomposable' means that there 
do not exist nonempty Ci and C2 s.t. the pairs Pi, P2 are always compatible for Pi G Ci, 
P2 6^2. The weight $c = Y^i-i^nv ciYlp^-p Pp^ 'is the sum over all monomials in the 
Taylor- expansion corresponding to multi-indices I with Ip > 1 for all P & C and c/ is 
the corresponding combinatorial factor, depending only on the incompatibility relation. 
Moreover, we have the decay- estimate of the form 

(A* \ 
yj <a*|P|, where a* 0.633 (2.16) 

for any fixed P. Here the sum is over all clusters incompatible with P, i.e. containing at 
least one polymer incompatible with P and we have put \C\ = X^pec l-^l "^here \P\ is the 
number of bonds of the polymer P. 



The proof is the same as that provided in [KOI]. Only the result was formulated for a 
translation-invariant setting, and applied as a technical tool in a different situation. (It 
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relies on the the general Kotecky-Preiss estimate [KP86] . A simpler proof of this kind of 
result is given in [BoZaOO].) 

We note that in our case A is estimated from above by 



A < A(a) := sup ^ gi {ux,y) < gi {u) ^^.n) 

where we have put 



u := sup Ux,y 
xeT 



(2.18) 

yeV-.y^x 



The second inequality of (2.17) follows from the positivity of the Taylor coefficients of 
gi. This estimate explains the occurence of the norm || ■ ||r. Such an estimate will be 
used over and over below. 

To compute the logarithm of the Laplace transform up to quadratic order in a we 
need only keep clusters with at most two bonds. We get from the general estimate on 
cluster sums provided by (2.16) the bound 



lo. 



g^(e^^("))-5:iog^(e^^)- ^in- E ^m,^.} 



xeFr- P:\P\=1,2 {^'l-P2}^|J'll = l^'2l = l.-Pl#-P2 

x{Pi)nx(P2)#0 

3 



<a*|r,| (^^j fori.<log(l + A*) 

(2.19) 

The hard part of the Taylor-expansion is now done by the general estimate. It remains 
to do some less elegant but elementary work: We still need to expand the three sums 
appearing under the modulus on the l.li.s. up to quadratic order in a, estimate the 
remainder terms and verify that they can be estimated in terms of the norms we have 
introduced. The quadratic order term obviously produces i^u {Xr{a))'^. 

Now, the first sum is trivially estimated. Let us define the function l{x) = — log(l — 
x) — X — X^fc=2 ^^/^- We have 



log/.(e^^)-l^(i;2) 
< l{92{v)) + gsiv) 

with 



= |iog(i + M^2(Ko)))-M^2(K.)) + M^3(K.))| ^2 20) 



V := SUp||ya;||oo f2 21) 

x&r ^ ' ' 



Here we used that /i(Va;) = 0. 

Let us come to the cluster sums. The cluster weights are obtained by comparing 
Taylor-coefficients (or by the inclusion-exclusion formula). One always has for single 
polymer clusters appearing in the second sum under the modulus of (2.19) that ^{p} = 
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log(l + pp). For these clusters we will write ^*{p} — (log(l + pp) — pp) + pp. Using 
l^{Ux,y) = we see that the activity of the single bond polymer P = {a;, y} is in fact 
of quadratic order in a. [This is better than the application of the bound (2.12) which 
holds for all polymers would show.] Indeed, 



\pp{ol)\ 



+ m{92{U.,y)) 



(2.22) 



< Ua;,ygi {2v) e^'" + g2 {ux,y) 



Therefore, log(l + pp) — pp is of fourth order, for both \P\ = 1,2. Thus we need to 
expand pp up to second order and control the third order error terms for both |P| = 1,2. 
Finally the ^s^p^^p^yierrcL is of forth order, too. To control its magnitude it is convenient 
to use again Proposition A.l using the improved bound (2.22). 

Now, let us give some more details on the estimation of the error terms. To estimate 
the difference between the cluster weights appearing under the second sum in (2.19) and 
the corresponding activities we use |${p} — Pp\ < 1{\pp\) to get 



P:|P| = 1,2 



< 



I {Ux,y9l i^v) + g2 {Ux,y)) 



+ J2 ^{9liUx,y)9liUy,z)) 



(2.23) 



Using the fact that all the functions appearing have positive Taylor coefficients we may 
estimate the r.h.s. by 



r,| Q/ {ugi {2v) e^" + (u)) +l(^gi{uf) 



(2.24) 



Next we need the error terms for the quadratic approximation on the polymer weights. 
Keeping the second order terms and using similar arguments as before we get for the 
single-bond polymer 



1 



pp{a) - p {Ux,y{Vx + Vy)) - -p (C/J J 



< Ux,y [2vgi {2v) + g2 {2v) e^"] + -ul^g^ {2v) (1 + e^^) + ^3 {ux,y) 



(2.25) 



For a double-bond polymer P = {{x,y}, {y, z}} we get in a similar fashion 



\pp{a) - p{Ux,yUy^z)\ < Ux,yg2{uy,z) + Uy^zg2{ux,y) + g2{ux,y)g2{uy,z) 
+ {ux,y + %,z) [^vgi {2v) + g2 {2v) e^""] + ]-Ux,yUy^zgi{'iv){l + e^"") 



(2.26) 
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Summing over the polymer and using the positivity of the Taylor coefficients of l,gi,g2 
we obtain 



^iog/.(e^^)+ Pipy-lf^iMc^f) 



P:\P\^l,2 



< \Tr\{l{g2{v))+g3{v) 



+ \u {2vg^ {2v) + ^2 (2^) e^"] + \u^g^ {2v) (1 + e^^) + \gz in) 
+ 1ug2{u) + g^iuf + 2« \lvg^ {2v) + g^ (2v) e^^] + \u^gx{Zv){\ + e^") 



(2.27) 



Let us finally treat the last cluster sum under the modulus on the l.h.s. of (2.19) involving 
two single-bond polymers. For a pair of incompatible polymers Pi, P2 one always has 
by the inclusion-exclusion formula that ^{Pi^pj} — log(l -|- pp^ -\- pp^) — log(l -|- pp^) — 
log(l + PP2). The easiest way to treat this term here is by application of Proposition 
A.l to the restricted polymer system that contains only single-bond polymers. We can 
use the improved second order bound (2.22). Denoting by the corresponding cluster 
weights we thus have from (2.16) that 



C-.CicpP 



|C| 



< a*|P| 



with the same a* ^ 0.633 and 



A' := sup V \p{x,y} I < ugi {2v) e^" + ^2 («) 

a;€r 



So we get 



E 



\^{Pl,P2}\ 



< 



a 



{Pl,J'2}:|Pl 1 = 1^21=1. ^'l 7^ -P2 
X(Pi)nX(P2)7^0 



(A*)2 



Collecting terms we arrive at the final estimate 

log/.(e^'^(«))-i^(X,(a)^) 



< \Tr\h{u,v) 



with 



h{u, v) 



a" 



(A*)3 



g^{uf + l{g2{v)) + gz{v) 



+ ^-u [2vg^ {2v) + g2 {2v) e^-] + ^-u'g^ {2v) (1 + e^-) + ^^3 («) 

+ 2ug2{u) + g2{uf + 2« [2^;^i (2^) + (72 (2^) e^"] + ^u'gi{?,v){l + e^") 



+ 



hi {ug^ {2v) e^^ + ^2 («)) + ^ (^i(«)') + ^ (2^) e^^ + ^2 («)) 



(2.28) 



(2.29) 



(2.30) 



(2.31) 



(2.32) 
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Now we use that u < 2K\\a\\randv < i^||a||r- So, K\\a\\r < |log(l+A*) =: d fa 0.05258 
implies that the cluster expansion is convergent. We get for these K||Q!||r the third order 
norm-estimate on the higher terms in the form 

h{u, v)<h (2K||a||r, i^||a||r) < D{K\\a\\rf (2.33) 

with D := s^\>x:0<x<d — ^'^^3''^'' ■ This is clear by the positivity of the Taylor- 

coefficients of h{x). It is then a trivial matter to compute the constant D given in the 
claim of the proposition numerically. We get -D <, ~ 4352 -|- 63 -|- 124 < 4540 where the 
first number gives a bound on the first term, the last number a bound on the last term, 
and the middle number a bound on the remaining terms of (2.32). This shows that the 
error term coming from the estimation of the higher order terms in the cluster expansion 
that depends on u alone provides by far the main contribution. 

Remark. Of course one cannot expect the series to converge without any smallness as- 
sumptions on ||a||r- In fact, for T = with v >2 and a{x) — J'\-\x\=i we are back to 
the usual ferromagnetic nearest-neighbor Ising-model, and the series is known to diverge 
for large J due to the existence of a phase transition. 



3. THERMAL DISLOCATIONS 

It is not too difficult to go through the proof given in the previous section to accom- 
modate the case of Model B of thermal dislocations. There are some changes, however. 
First of all, we need a different norm estimating the variation of the (Fourier transform 
of) the observable w.r.t. variations up to the magnitude 8. We define the semi-norm 



|Q^||r,(5 := sup y sup \a{x — y + z) — a{x — y + z')\ 

1% \^\Az'\<25 



(3.1) 



Note that this seminorm vanishes on constant functions. We denote 

Yr{a) = \Tr\ {i^{a) - {^^{a))) (3.2) 

for the nonnormalized centered autocorrelation measure applied to the function a. 

Then we have a norm-estimate on the Laplace-transform that is analogous to Propo- 
sition 1. The result is the following 

Proposition 2. Suppose that the dislocations u) = {wx)xer have independent, not nec- 
essarily identical distribution /i, such that \ujx\ < S < 00, for all x fj,-a.s. 
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Then there are universal constants d > 0, and D < oo, independent of the set T and 
the distribution ^, such that, whenever a :W ^ 'C is small enough such that ||Q:||r,5 < d 
and \oi{x)\ = \a{—x)\ for all x, we have the estimate 



< \Tr\D \\a\\^^g 



(3.3) 



The values of the constants can be chosen like d = 0.0525 (same as in Proposition 1), 
and L> = 4.38 •103(< D). 

Proof. We follow the lines of the proof given in the previous section. Our 'Hamiltonian' 
now becomes 



{x,2}cr^ 



(3.4) 



with the pair potential 



-\-a{y - X -\- ujy -uj^ - ii[a{y - X -\- (jjy -a;^;)^ 



(3.5) 



for X ^ y. Note that there is no single site potential this time, since the corresponding 
expression vanishes for a; = y. We note that 



\U. 



x,y\\oo 



<2 sup \a{x — y -\- z) — a{x — y -\- z')\ 



z,z' et 



\z\,\z'\<15 



(3.6) 



So we have 



u' := sup ^ Ux,y < 2||Q!||r,5 ^3 7) 

Now the steps of the proof of Proposition 1 stay true, leading to formula (2.31) with 
v = 0. 



< \rr\h{u',v = 0) 



(3.8) 



with the function h given in (2.32). The constant d stays the same and for the constant 
D we get the better value D = ^^^f"^ <^ 4352 + 10 + 12 < 4380. This shows that we 
get essentially the same constant as that of Proposition 1 and the diagonal terms didn't 
do much harm. 
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4. MORE BOUNDS, CLT, AND FINAL PROOFS 

In this chapter we state the more detailed versions of Theorems 1 and 2, along with 
their proofs and also provide the Central Limit Theorem. In particular the results will 
still contain the discrete F-dependent norms (2.2) rasp. (3.1). We have preferred here to 
write the estimates in terms of the autocorrelation-measure applied to a function (rather 
than its Fourier-transform) . We will use here the notations from Chapter 2 for Model A 
and Chapter 3 for Model B. 

Let us start with the result for Model A. Suppose that ct is a given function on M.'^ with 
|Q!(a;)| = |q;(— a;)| for all x e W^. Recall the definition of the discrete norm ||Q;||r given 
in (2.2). Recall the notation Xr{a) = \Tr\ (7^(ct) — //(7^(a))) for the nonnormalized 
centered autocorrelation measure applied to the function a. In this situation the following 
result holds. 

Addition to Theorem 1. We have the large deviation estimate 



fx(\Xr{a)\>s\Tr\) <2expl -|F^| X ^;s,n (4.1) 



for all e>0, where Sr = (^X^ {k^^)) ^ 4. 

Here, for fixed s > the function jd,D{ ■;•§)• [0; oo) ~^ [0; oo) has the form 

■ ^ \ ^{{l2De + s'')i-lSDes-sA, ife<d{s + ZDd) 

Md(£;s) = < \ ^ (4.2) 

I. d {e - d ^- dD)) , else 

where d > 0, D < oo are the same numerical constants as in Theorem 1. It is convex, 
nonnegative, strictly increasing in e. It is decreasing in s and in D, and increasing in d. 

Remark. The statement is stronger than the simpler one given in Theorem 1 in two ways. 
First of all, the Sobolev-norm || ■ \\jy^a appearing therein is replaced by the sharper norm 
II ■ ||r introduced in (2.2). This is only very minor because the Sobolev-norm will be used 
in practical applications. Next, we have kept the normalized variance s^- Usually will 
be of the order unity, e.g. when F is a lattice and the scatterers are i.i.d. There can 
however be cases of sets F^. and functions a for which this quantity will go to zero with 
r ] oo. 

Now, putting together the two pieces of information < 4 and || ■ ||r < || ■ \\v,a (see 
Appendix B, Proposition 3) the simplified statement of Theorem 1 immediately follows, 
by the monotonicity of the function j in s. 

The situation is completely analogous for Model B of thermal dislocations. To for- 
mulate the corresponding statement we recall definition (3.2) for the non-normalized 
autocorrelation-measure applied to a. Recall the discrete F-dependent semi-norm (3.1). 
Then our result reads as follows. 
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Addition to Theorem 2. Suppose that S < oo. Then we have 

l^(\Yr{a)\>e\Tr\) <2exp||-|r,| x j^^^ (p^; g,) j (4.3) 

for all e > 0, where qr = ji^A* (^r ^^^^'^ d > 0, D < oo are the same 

numerical constants as in Theorem 2 and the function fy is given in (4-2). 

Remark. The statement of Theorem 2 follows from here by g^ < 4 and the norm estimate 
given in Appendix B, Proposition 4. 



Looking at the variable on the central-limit scale we get the following result. 

Theorem 3. Suppose that lim^i-oo A* (^r ('^)) l^rl"^ — oo. Then the standardized vari- 

_ 1 

able Xr{a) [X^{a))) ^ converges weakly to a standard Gaussian distribution. The 
same statement holds for replacing X^.. 



Finally we give the proofs. 

Proof of the Addition to Theorem 1. Assuming the uniform estimates on the Laplace 
transform provided in Proposition 1 it is a trivial matter to derive the large deviation 
upper bound. Indeed, by the exponential Chebychev inequality we have 

^i{±Xr{a) > \rr\e) < inf e'^fel^-^'V (e^^'^^*"/^^""!!^))) 

t:0<t<d V / 

<exp('-|r,|x sup f--^t--^fi(x^(-^))i--Dt^^^ 
\ t:0<t<d\K\\a\\r \Tr\ \ \K\\a\\rJJ 2 



Call the function appearing in the exponent in the bound 

jd,D(£; s) := sup (et-^- DtA (4.5) 

t:0<t<d V ^ / 

Observe that jd/\,x^D{^^', ^"^s) = jd,D{£'iS). A simple computation shows then that j 
has in fact the explicit form (4.2) given in the Addition to the Theorem 1. (In the small 
£-range the maximizer is t = -^+^^^^+^^ ^ the large e-range the maximizer is t = d.) 
j is convex and nonnegative. The monotonicity properties claimed in the Addition to 
Theorem 1 are now immediate by formula (4.5). Finally, it is a simple exercise to see 
that Sr < 4, using the independence of the scatterers. 

The proof of the Addition to Theorem 2 is the same. 
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Proof of Theorem 3. Prom the bound on the error for the quadratic approximation of the 
Laplace transform given in Proposition 1 follows immediately that, for all fixed t e C we 

have that limj.-foo log A* (^^xp (^Xr{a) [fj, [X^{a))) ^ = ^^/2, under the assumptions of 
the theorem. This shows the claim. The proof for Model B is identical. 



APPENDIX A: SCATTERING THEORY FOR POINT SCATTERERS 

Let us briefly recall the basic elementary formulae of scattering theory that describe 

the connection between the autocorrelation 7^ (resp. 7^ ) and the scattering image. For 

notational concreteness we only consider Model A (fixed locations). (For more on this see 

[Hof95b], [BaaHoeOO]). Suppose a beam with wavelength A;, hits the finite collection of 

point-scatterers located in the finitely many points F^. Denote by cq e M*^ the incoming 

direction (where |eo| = 1 is a unit vector). The modulus of the scattering amplitude rjx 

gives the amplitude of the scattered wave and the phase of r]x gives a local phase shift at 

the site x. Consequently the intensity of radiation scattered elastically in the direction 

2 

e is given by X^xer^ ^e*'^'^ with k = 27r(e — eo)/Ab. To understand the l.h.s. of this 
formula take ry^; = 1 and note that in this case k ■ x is the phase difference of a beam 
scattered at site x relative to that of a beam scattered at a (hypothetical) scatterer at 
site 0. 

Multiplying the intensity by a test-function (f{k) (that models the sensitivity of a 
counter) and normalizing by the number of scatterers then leads to the quantity 7^ {(f) = 
7^ (</?). Here we choose the convention (p{x) — J^^ e^^''^(p{k)dx to define the Fourier- 
transform of a Schwartz- function (p. The Fourier-transform of a tempered distribution is 
then defined by duality. 

So, when one is interested in the infinite volume limit, one likes to look ([Hof95b]) 
at the scattering measures 7^ in the sense of (tempered) distributions and is interested 
in the weak limit r | 00, i.e. lim^ioo 7r (V') where </? is a Schwartz function. Then, if a 
limiting distribution exists at all, it can have a discrete part, an absolutely continuous 
and a singular continuous part, the discrete part (Bragg peaks) caused by 'order', the 
continuous parts showing diffuse scattering caused by 'disorder' of the scatterers. This 
is in analogy to statistical mechanics where sharp phase transitions occur only in the 
infinite volume. 

Disorder-averages of diffraction measures 

Our theorems give us good control over the fluctuations of the scattering measures 
7^. The estimates are independent of the behavior of the mean, and the nature of the 
limimiting distribution, if it exists. To compute their disorder averages of the scattering 
measures is a trivial matter. We get 
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Model A: J fi{dvW{^) = + ^ E.er. (^(I^-I') " Hv.W) x ^(0) 

where n{'r]x) —'■ mx is the mean-vahie of the scattering amphtude and m = {'mx)xer- The 
first term describes the scattering image of a system where the scattering amphtudes have 
been replaced by their means and the second term a homogenous diffuse background. 
So we see, that a.s. convergence for the averaged scattering images holds if and only 
if the two individual terms converge. This is true for F a crystal or quasicrystal and 
rjx are i.i.d., with e.g. Fr being increasing balls. The latter statement follows since a 
crystal or quasicrystal is known to possess a natural autocorrelation function. To see 
how to construct an example of independent but not identically distributed scatterers 
on a quasicrystal for which the mean converges, see Paragraph 7 of [BaaMoo98]. On 
the other hand, it is simple to construct examples of systems on lattices with prescribed 
convergence/non-convergence of each of the two terms along a given sequence of volumes. 
This is done by choosing the distribution of ?7a;'s in a non- homogenous way; think of a 
sparse sequence of increasing volumes Tr and choose two different distributions in the 
annuli rr+i\rj. for r even resp. r odd. Still, also in these examples without convergence of 
the mean, under the assumption of uniform boundedness of the distribution, selfaveraging 
in the sense of Theorem 1 would hold. 

Model B: / fx{duj)^l;! {(p) = ^ E^^^'er, (^(^ - x' + ujx - ujx')) + ^(0) 

Again, one can construct artificial distributions of dislocations of scatterers on a lattice 
such that this expression does not have a well-defined limit. Choose e.g. a;^ = for x 
in the annuli rr+i\rr for r even, and a non-trivial bounded law for ujx for r odd. Still, 
self-averaging holds. 

However, if the cj^j's are i.i.d. with single-site distribution we get from this 
/ n{duj)^!^{k) = 7°(A;)|/i(A;)p + (l — where 7°(/c) is the density of the Fourier- 

transform of the autocorrelation with all the scatterers sitting at their sites in Tr and 
= |//(e*'^'^ '^)p is the famous Debye-Waller factor reducing the intensity of the 

reflexes. 



APPENDIX B: NORM ESTIMATES 

Finally we like to give the norm-estimates along with their proofs that are needed to 
obtain the final form of the theorems as they are stated in the introduction. 

Proposition 3. Suppose that FcM'^ and the number a is a bound on the minimal dis- 
tance between the points in F. Then we have the bound on the discrete T-norm in terms 
of the a-weighted Sobolev-type norm of the form \\g\\r < IblU.a- 

Remark. By scaling one may construct examples that show one can not do with less than 
the first 1/ derivatives, in general. 
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Proof. Put disjoint balls of radius a/2 around the points of F. Consider anyone of them. 
And assume without loss that its center is 2; = 0. Write for simplicity B = 5^/2 (0). 
Then we have |^(0)| x \B\ < Jg \g{y) ^ g{0)\dy + Jb To express the first integral 

on the r.h.s. as an integral of the derivatives of g over B we use Polar coordinates 

use the one-dimensional Taylor-expansion of the 
function r 1— g{re) —: Xe{r) of the radial coordinate r up to order u — 1. Expanding 
around the point r we get \g{re) - g{0)\ < Y^lZl i\Xe''\r) \ + dsjf^\xP {s)\. This 

gives /g \g{y) - g{0)\dy < Y,k=o ''"If / ^{de) Jq^^ dr r''~'^\xi^\r)\. The reader should 
check that also the term ior k = u can be bounded in this form (interchange the orders of 
integration between s and r!) This argument only works since the power v — 1 reappears 
under the integral of the remainder term. Dividing this inequality by the volume of B, 
bounding the A;-th directional derivatives by and integrating over the whole of 

now proves the claim. 



For the semi-norm ||(7||r,5 that was introduced in (3.1) (needed to control Model B) 
we get the following analogous estimate. 

Proposition 4. Suppose that FcM.^ and the number a is a bound on the minimal dis- 
tance between the points in F. Assume that a := a — AS > 0. Then we have the bound in 
terms of the a-weighted Sobolev-type semi-norm 

\\9\\r,s < ELo i.wrn^ I^'^XB.MO) \\d'^'9{y)\\dy (< m9L,a 



■g{x-y + z')\ < 
g{w + te) . 



Proof. For fixed x ^ y vnT and any \z\, \z'\ < 25 we have \g{x — y + z) 
!o sup|e|=i I f^g{x -y + z^te)\dt<A5 sup^^e^iCa^-y) sup|e|=i it 

Using the estimate in terms of the integrals over balls in terms of derivatives up to 
order of the dimension provided in the proof of Proposition 3 we get for w e B2s{x — y) 



t=o 



g{u + te)||(itt 



3ere we have used the radius a/2 because this implies that Ba/2{w)cBa/2{x — y), 
independently of w, and so we get that the r.h.s is bounded by 
j^Tl Efc=o hja/W^ \\d''+'^g{u)\\du. This gives the desired estimate by sum- 

ming over y that are not equal to x, and extending the integral over all of M.^\Ba/2{0)- 
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